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Abstract 

We study the creation of photons in a one dimensional oscillating cavity 
with two perfectly conducting moving walls. By means of a conformal trans- 
formation we derive a set of generalized Moore's equations whose solution 
contains the whole information of the radiation field within the cavity. For 
the case of resonant oscillations we solve these equations using a renormal- 
ization group procedure that appropriately deals with the secular behaviour 
present in a naive perturbative approach. We study the time evolution of the 
energy density profile and of the number of created photons inside the cavity. 

I. INTRODUCTION 

It is well-known that in the presence of moving boundaries the vacuum state of the 
electromagnetic field may not be stable, which results in the generation of real photons. 
The generated radiation exerts pressure on the moving boundaries which can be looked 
upon as a dissipative force that opposes itself to the mechanical motion of the boundaries. 
The generation of photons, which is an amazing demonstration of the existence of quantum 
vacuum fluctuations of QED, is referred to in the literature as dynamical Casimir effect |]| 
or motion-induced radiation 0. It goes without saying that it would be very nice to have 
an experimental verification of this prediction. Due to the technical difficulties involved in 
the detection of the phenomenon, up to now no concrete experiment has been carried out, 
and there are only some few experimental proposals ||(§]. However, a feasible experimental 
evidence is not out of reach, and therefore it is of interest to explore different theoretical 
models to describe the process and identify signatures which permit to distinguish vacuum 
radiation from spurious effects. 

Research in the field has mainly concentrated on one dimensional models, which are 
useful for giving an account of the main physical processes participating in the phenomenon 
(a small number of works deal with more realistic three dimensional models 0-0). In 
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this work we will also restrict ourselves to one dimesional models. Motion induced effects of 
vacuum radiation already show up for a single mirror moving with a non-uniform acceleration 
in vacuum |7|]. Since the amount of radiation generated is very small, basically determined 
by the ratio of the speed of the mirror to the speed of light, much attention has been paid 
to the study of one dimensional models for which the effect is enhanced. 

A cavity made of two perfectly parallel reflecting mirrors, one of which is motionless and 
the other one oscillates with a mechanical frequency equal to a multiple of the fundamental 
optical resonance frequency of the static cavity is a thoroughly studied example where such 
an enhancement takes place [JS-H]. It is typically considered that the cavity is motionless 
and that at some instant one mirror starts to oscillate resonantly with a tiny amplitude. For 
small times after the motion starts, one can make a perturbative expansion of the equations 
of motion of the field in terms of the small amplitude to find an approximate solution. In 
this way one can study the structure of the electromagnetic field inside the cavity, which 
departs from the standard static Casimir profile (which is constant over the whole cavity) 
and develops a structure of small and broad pulses. The number of mot ion- induced photons 
grows quadratically in time, and the spectrum has an inverted parabolic shape with an 
upper frequency cut-off given by the mechanical frequency, its maximum being at half that 
value PUTS]. Similar results are found in || by means of a scattering approach for the 
radiation emitted out of a lossy cavity. However, for long times these methods are not valid, 
and new approximation techniques are required. In |13],[16| it is shown that in such a limit 
the structure of the electromagnetic field is non trivial, with a number of pulses equal to 
the mechanical resonant frequency, whose width decreases exponentially and whose height 
increases exponentially with time, in such a way that the total energy within the cavity 
grows exponentially at the expense of the energy given to the system to keep the mirror 
moving. Also, the spectrum does not have an upper frequency cut-off. Through a process of 
frequency up-conversion, the generated photons contain frequencies of higher order cavity 
modes and thus exceed the mechanical frequency. The physical mechanism of such an optical 
pumping into the high frequency region is the Doppler up-shift of the field upon reflection 



at the mirrors. Similar conclusions are found for the lossy cavity (17| . 

The case of cavitites with two moving mirrors has also been considered recently. In 
the small time approximation, both for the ideal cavity |Tj| and for the lossy one 0, it is 
found that the number of motion-induced photons grows quadratically in time and that the 
spectrum is once again parabolic. In the long time approximation, the lossy cavity has been 
studied with the scattering approach jjI7 |. Just as in the case of a single oscillating mirror, 



it is found that in this regime, and for the two types of motion described above, there is 
pulse shaping in the time domain and frequency up-conversion in the spectrum of emitted 
photons from the cavity. A striking feature of the spectrum is that no photons are emitted 
at frequencies equal to multiple integers of the mechanical frequency. 

In this work we will consider an ideal cavity with two mirrors oscillating resonantly at 
the same frequency, and we will allow for different amplitudes and a possible dephasing 
between the mirrors. To investigate the problem we will deduce a generalization of Moore's 
equation || to the problem of two moving boundaries, whose solution gives the complete 
information on the electromagnetic field inside the cavity. For the case of resonant harmonic 
motions, no exact solution exists, and we find an approximate analytic solution based on 
a renormalization group (RG) technique. We have already applied this method in [IB| for 
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the case of a single oscillating mirror, and just as in that case, the strategy allows to find a 
single solution valid both for small and long times. This will allow us to describe precisely 
the bahaviour of the energy density and the number of photons for all times. 

As we shall see, motion-induced radiation strongly depends on the relation among the 
amplitudes of oscillation, the frequency and the dephasing. For some relations among these 
variables, there is constructive interference and a series of pulses develop within the cavity 
that grow exponentially in time, and frequency up-conversion takes place. For some other 
relations, there is destructive interference and hence no vacuum radiation. We also show 
that our solution is capable of accounting for other physical behaviours, for which the peaks 
grow quadratically rather than exponentially. 

The paper is organized as follows. In Section 2 we will introduce the generalization of 
Moore's equations for a moving cavity and we will explain how to calculate the energy density 
and the number of motion-induced photons. In Section 3 the renormalization group method 
is described and applied to the problem of harmonically oscillating walls with dephasing. 
In Section 4 we study some particular dephasings, which we will call translational and 
breathing modes. In Section 5 another motion is considered, which has a cualitatively 
different behaviour as compared to those of the previous Section. Finally in Section 6 we 
make our conclusions. 



We consider a one-dimensional cavity formed by two perfectly reflecting mirrors, each 
of which follows a given trajectory, say Lit) for the left mirror and R(t) for the right one. 
These two trajectories are predetermined (i.e. are given data for the problem) and act as 
time-dependent boundary conditions for the electromagnetic field inside the cavity. The field 
equation for the vector potential takes the form of the equation for a massless scalar field 
(— dt +dl)A(x, t) = Q and the boundary conditions are A(x = L(t), t) = A(x = R(t),t) = 
for all times. If we express the field in terms of creation a k and annihilation operators 
for photons in the form 



then the mode functions tpk{x,t) must be chosen so as to satisfy the above boundary condi- 
tions. 

In the case when only one of the walls moves, say the right one, the modes can be written 
in terms of a function U(t) as 



II. GENERALIZED MOORE'S EQUATIONS 
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(1) 
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(2) 



The speed of light is set to unity. 
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and the boundary condition on the right is met provided that the function U verifies 
U(t + R(t)) — U{t — R(t)) = 2, which is known as Moore's equation ||. The complete 
solution to the problem involves finding a solution U (t) in terms of the prescribed motion 
R(t). Moore's equation can in fact be deduced by means of a conformal transformation from 
the original spacetime coordinates (t, x) to a new set of coordinates (t, x) in which not only 
the left mirror but also the right one is fixed. This transformation takes the form 

i+ x = U(t + x) , i-x = U(t-x), (3) 

which, after mapping the coordinate of the left mirror as L(t) = — > xl = and for the 
right one R(t) — > xr — 1, leads to Moore's equation. 

Let us now consider the more general case in which both mirrors move. Evidently we can 
also make a similar conformal transformation, but know we need two functions U instead of 
one. Defining the transformation as 

t + x = G(t + x) , t-x = F(t-x), (4) 

and mapping L(t) and R(t) as before, we obtain a set of generalized Moore's equations 

G(t + L(t)) - F(t - L(t)) = 
G(t + R(t))-F(t-R(t)) = 2, (5) 

which, when solved for given motions for the mirrors, allow us to find the solution for the 
modes inside the cavity. Indeed, the modes can be casted in the form 

t) = [ e -^ G ^ - e -^F(t- x) i ) (6) 

and they satisfy both the field equation and the boundary conditions. 

We shall be interested in studying the spacetime profile of the energy density of the field 
between the moving walls 



(T 00 (x,t)) = 




where the expectation values are taken with respect to the vacuum state. As is well-known, 
this quantity is divergent and a regularization method is needed to get meaningful results. 
Using the point-splitting method and introducing advanced u = t + x and retarded v = t — x 
coordinates, the energy density can be rewritten in terms of the functions G and F in the 
following way 

oo 

(T 00 (U,V)) = ~Y, k {G'(U)G'(U + ie ) e -^lG(u)-G(u+ie)} + 
k=l 

F'(v)F'(v + ie ) e -ik*[F(v)-F(v+U)] J ) (g) 



2 The boundary condition on the left (fixed) mirror is automatically fulfilled by this form for the 
modes. 
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with e — > + . From here it is straightforward to get the renormalized version, (Tqq(x, t))ren 
-f G {t + x) - f F (t - x), where 
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Suppose that for t < both walls were at rest separated by a distance A and that the 
field was in its vacuum state. The solution of the generalized Moore's equations is simply 
G(t) = F(t) = t/A and the mode functions ipk correspond to positive frequency modes. If at 
t = the boundaries begin to move, it is well-known that for some types of motion the field 
does not remain in vacuum, but photons are produced through nonabiabatic processes. A 
consistent calculation of the number of created photons through mot ion- induced radiation 
requires to have a well-defined vacuum state in the future. To this end we consider that 
at time t = T both walls come to rest (L(t) = and R(t) = A for t > T). The evolved 
old vacuum state does not coincide with the vacuum in the future, but rather it contains 
a number of photons, which can be calculated by means of the Bogoliubov coefficients 
Arm = (^n>^™ )*> where ijjm\x,t) = (7rm) -1 / 2 sin(m7rx/A)e~ i7rm * //A is the mode function 
for the static problem and ip n (x,t) is the mode function which solves the nonstationary 
problem for t > and coincides with ipn\x, t) for t < 0. Writing the mode functions in 
terms of the functions G and F, integrating by parts, and using the set of Moore's equations 
to drop the surface terms, we get the following relation for the Bogoliubov coefficient 



1 fm \ f t/A , r . 

— \ — < / dxexp{—z7i[nF{Ax) + mx\\+ 

2 V n I Jt/A-l 



t/A+l 



dxexp{— iir{nG(Ax) + mx]} 



(10) 



t/A 



for times t > T. The number of created photons inside the cavity after the stopping time T 
in the n-th mode is given by N n (T) = Ylim \Pnm(t,T)\ 2 (the dependence of the Bogoliubov 
coefficient on t is just a phase) and summing over n we get the total amount of motion- 
induced photons. 



III. RENORMALIZATION GROUP METHOD FOR MOORE'S EQUATIONS 

For resonantly harmonic motions it is not possible to find an exact solution to Eqs.(0) 
and approximation methods are compelling. A naive approach is to make perturbations 
in the amplitude of the oscillation, but it turns out that the strategy is ill-fated, because 



3 The inner product is the 



usual for Klein-Gordon equation, namely 
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of the appearance of secular terms proportional to the time which, after a short period, 
make the approximation break down. In |16| we have applied a method inspired in the 



renormalization group to treat these singular perturbations for the case of a one dimesional 
cavity with one oscillating mirror. The method has a wide range of application in different 
fields, especially for studying ordinary differential equation problems involving boundary 
layers, multiple scales, etc. fl9j . In the following we shall extend the method for the case 
of a one dimensional cavity whose mirrors oscillate in resonance with the cavity. More 
specifically, we consider that for t < the two mirrors are motionless and separated by a 
distance A, and that at t — they start to move as 

L{t) = eA L sin = + e5L(t) 

R(t) = A - eA R sin (0) + eA R sin + <pj = A + eSR(t), (11) 

for the left and right mirrors respectively. Here is a possible dephasing angle, A^ and A R 
are amplitudes of oscillation and e C 1 is a small parameter. 

Let us first start with the perturbative approximation. We expand both unknown func- 
tions G(t) and F(t) in terms of the small parameter e and retain first order terms only, 
G(t) = Go(t) + eG\{t) and F(t) = F (t) + eFi(t). Equating terms of the same order in the 
set of generalized Moore's equations, we get for the zeroth order part 

G (t) - F {t) = (12) 
G (t + A)-F (t-A) = 2, (13) 

and for the first order part 

G l (t) - F x (t) = -6(t)5L(t) [G' (t) + F '(t)] (14) 
Gt(t + A) - Fx(t - A) = -6(t)8R(t) [G' (t + A) + F^(t - A)] . (15) 



The general solution to Eqs.fll^, |13|) is 



G (t)=F (t) = c+- + }_^ Ansml — J +B n cosf — J 

n>l L V / V / 



(16) 



where c, A n , and B n are constants determined by the boundary conditions of the problem. 
These are obtained from the fact that the modes ^(x, i) must be positive frequency modes 
for t < 0, which implies that G(t) = t/A for -oo < t < A and F(t) = t/A for -oo < t < 0. 
The different ranges for the functions G and F follow directly from the nonlocal structure 
of Moore's equations. 

Making the shift t — > t — A in Eq.flT4|) and replacing the result in Eq.(^) we get an 
equation for the first order correction to the function G, namely 

G x {t + A) - G x {t - A) = 6{t - A)8L{t - A) [G' (t - A) + F^t - A)] - 

e(t)5R(t)[G' (t + A) + F^t-A)]. (17) 

Since this equation is linear, the solution is of the form G\ = G± + G± , where 
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G { i\t + A) - G^\t - A) = 6(t - A)5L(t - A) [G' (t - A) + F '(t - A)] (18) 
GS 2) (t + A) - GS 2) (t - A) = -6(t)5R(t) [G' (t + A) + F^t - A)] , (19) 
whose general solutions read 

G?{t) = ^ t -6{t) sin ifpct/ A) x 

1 + vr ^ n [A n cos (n-Kt/A) - 5 n sin (rort/A)] I + # (1) (t), (20) 

n>l J 



for Eq.(0), and for Eq.(|i~9"|) we get 

Ant 
A A 



G?\t) = -JL-9(t + A) {sin(0) + (-l) 9+1 sin(g7rt/A + 0)} x 



1 + tt ^ n[A n cos {mrt/ A) - B n sin (nvrt/A)] I + # (2) (t), (21) 



n>l 



where and g( 2 ) are arbitrary periodic functions of period 2A. The first order correction 
for F can be deduced from the first order correction for G we have just found using Eqs.(|T4"D 



or (15) indistinctly. We see that the perturbative corrections contain secular terms that 
grow linearly in time. Therefore the approximation will be valid only for short times, that 
is, et/A < 1. 

In order to determine the two unknown periodic functions we have to consider the bound- 
ary conditions for the functions G and F. We have already said that the nonlocal structure 
of Moore's equations implies that, although at t — the motion of the mirrors starts, the 
expression for G for times up to t = A is given by the solution for motionless walls, namely 
G(t) = t/A for t < A, while that for F reads F(t) = t/A for t < 0. If we assume that these 
boundary conditions are already satisfied by the zeroth-order solutions Go(t) and Fo(t), then 
the periodic functions must be chosen so that G\(t) = and Fi(t) = in the respective 
intervals. This fact, 
boundary conditions 



intervals. This fact, when translated to the functions and G^ imply the following 



G { i\t) = 0,for < t < 2A 

Gf\t) = 0, for -A <t < A, (22) 
which leads to the following expressions for the periodic functions 

g {l \{2k + 1)A + z) = -^^^ sin (qnz/A) x 

1 + 7r n [A n cos {nixz/ A) — B n sin (nnz/A)] > , (23) 



n>l 



and 



g {2) (2pA + uj) = -^^{sin(0) + (-l) 9+1 sin(g7rcu/A + 0)} x 



1 + 7r V^n [A n cos (niruj/A) - B n sin (jvku/A)] > , (24) 



n>l 
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where t = (2k + 1)A + z, k = 0, 1, 2, . . . and —A < z < A for the function g^ 1 ', while for the 
function we have t = 2pA + u, p = 0, 1, 2, . . . and —A < uj < A. Given t, the values of the 
integers k and p are obtained as k = p = |int(t/A) for hit (t/A) even, and k = |[int(t/A)-l], 
p = ~[int(i/A) + 1] for int(t/A) odd. Note that during the first period k = (p = 0), the 

function (g^) makes G± (G^) vanish identically. As we have already seen, since the 
mirrors were at rest for t < 0, we must impose G(t) = t/A for t < A and F(t) = t/A for 
t < 0. Therefore c = A n = B n = 0, and the perturbative solution for t > 0, to order C(e 2 ), 
is 



G(t) = j + sin (qnt/A) + [sin (0) + sin (qnt/A + 0)] (25) 

F(t) = G(t) + 2e^ sin (g?rt/A) . (26) 

These perturbative solutions suffer from the aforementioned secularity problems, being valid 
for times t/A <C e _1 . In order to deal with this drawback and get improved solutions valid 
for longer times, in the rest of this Section we describe the RG method we mentioned before 
for this problem of an oscillating cavity. 

What the renormalization group method does is to improve the perturbative expansion by 
resumming an infinite number of secular terms. In general, if one performs the perturbative 
expansion to higher orders, there appear different time scales, et to first order, e 2 t and e 2 t 2 to 
second order, e 3 t, e 3 t 2 and e 3 t 3 to third order, and so on. The RG technique sums the most 
secular terms of each order (e n t n ), and it is therefore valid for times t/A <C e~ 2 . The way to 



carry out the resummation is nicely described in |T9[ and it basically consists in introducing 
and arbitrary time r, split the time in the secular terms of the first order perturbative 
corrections as t = (t — r) + r, and absorb the terms proportional to r into the "bare" 
parameters of the zeroth order perturbative solution, thereby becoming "renormalized" . 
Introducing the arbitrary time r and splitting t as said, the perturbative solution can be 
written as 

t - T 



G(t, t) = c(r) + V" [A n (r) sin(mrt/A) + B n {r) cos(mrt/A)] + 



n>l 



A 



+e - T A ^ sin (qnt/A) I 1 + vr ^ n [A n (r) cos (nnt/A) - B n (r) sin {nut /A)} 



n>l 



t-T A 



+ el T^ ( sin M + sin (^/ A + 0} 



A A 



x 



1 + 7r n [A n (r) cos (nnt/A) - B n (r) sin (rwrt/A)] 



n>l 



+eg W(t,T)+egW(t,T), (27) 

where the bare parameters c, A n and B n have been replaced by their renormalized counter- 
parts c(r), A n (r) and B n (r). Here g^(t,r) and g^ 2 \t,r) respectively denote the functions 
<7^(£) and g^ 2 \t) with the same replacement. Note that these functions are no longer peri- 
odic due to the RG-improvement. 

Since the time r is arbitrary, the solution for G should not depend on it, which implies the 
following RG equation (dG/dr) t = 0. In our case it consists of three independent equations 



S 



<9c(r) 

dr 
dA n {r) 

dr 



dB n (-, 
dr 



\ + -{-l) q+1 b + 0^) 

A 7T 

-aS nq -2(-iy +1 bnB n 

71 

+ \n - q\ [aA\ n _ q \ - 6sg(n-q) J B| n _ f/ |] 
-(n + q) [aA n+q + bB n+q ] + 0(e 2 ) 
2 



(28) 



(29) 



-65, 



7T 



nq 



2(-l) q+1 bnA r , 



+ \n - q\ [asg(n - q)B\ n „ q \ + bA\ n - q \ 
+(n + g) [-aS n+9 + 6A n+g ] + C(e 2 ) 



where 



a 



I, 



e 7i 
A2 

A2 A 



^ + ^(-l)^cos(0) 



sin(0). 



(30) 

(31) 
(32) 



These parameters a and 6 play a crucial role because they determine the behaviour of 
the solutions to the set of generalized Moore's equations. There are four distinct cases. The 
simplest one is for a = b = 0, which happens, for example, for equal amplitudes A L = A R , 
zero dephasing and even frequencies. In this case there is no secular behaviour at the level of 
the perturbative solutions Eqs. (|25|j26D , which are then valid also for long times. The energy 
inside the cavity oscillates around the static Casimir value and there is no motion-induced 
radiation. A second case is a ^ and 6 = 0, which occurs, for example, for equal amplitudes, 
zero dephasing and odd frequencies. In this case secular terms do appear in the perturbative 
solutions and the RG method is useful for finding the long time behaviour, which shows an 
exponential increase of the energy in the cavity and motion-induced photons. This case will 
be the subject matter of the next Section. A third case is a = and 6 ^ 0, which takes 
place, for example, for a static left mirror Al = and dephasing = tt/2. Here there are 
also secular terms at the perturbative level, and for long times the energy does not grow 
exponentially but quadratically, photons also being generated. We shall deal with this case 
in the Section V. Finally, the case a ^ and 6 ^ is similar to the second case, there is 
motion-induced radiation and an exponential increase of the energy. We shall not cover this 
case in detail since the expressions for the solutions to Moore's equations are cumbersome. 



Now we solve the RG equations (p8| 



The solution for c is trivial, c(r) 



+ 



A A 5111 
Bn + B- 



t + n, with k a constant to be determined. Writing A r 
the new variables satisfy 



A n - A^ n and B n 



dAJj) 
dr 



2 

7i 

+ (m - q) 
— {m + q) 



ad mq - 2(-l) q+1 bmB r , 



dA m — q bB m _ q 



+ bB, 



dB m (r) 
dr 



-b5 ma + 2(-l) q+1 bmA 



m+q 



+ 0(e 2 



(33) 



7T 



J mq 
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+(m-q) aB m - q + bA m „ q 

+ (m + q) -aB m+q + bA m+q + C(e 2 ). (34) 

The initial conditions for these differential equations are dictated by the perturbative 
solution: c(0) = A m (0) = B m (0) = 0. This implies that the constant k is zero. In order to 
solve for A m and B m , we first decouple the equations through the transformation C m = A m — 
iB m and D m = A m + iB m , and introduce a generating functional M(s,t) = ^2 m s m C m (T). 
It is easy to see that this functional verifies the following differential equation 

?M. = -U - t b)s q + \-2i(-l) q+1 bs + (a - ib)s q+1 - (a + ^s 1 " 9 ] (35) 
or 7r L J os 

with boundary condition M(s, r = 0) = 0. The solution can be obtained by proposing 
an Ansatz M(s,r) = §[e~ T a(s)] + (3(s), where $[•••], oc(s) and @(s) are functions to be 
determined. We shall not dwell on the details of finding these functions, but suffice it to 
say that the last two are straightforwardly derived after introducing the Ansatz in Eq. (|35|) , 
while the first one is got once the initial condition on M is imposed. The solution reads 

2 

M(s, t) = - {i(~l) q+1 qbr - In cosh(gar)- 

7T 

In [(1 + i{-l) q+1 {b/a) tanh(gar)) - (1 - ib/a) tanh(gar)s 9 ] } . (36) 

Expanding this solution in powers of s (and doing the same for its complex conjugate) we get 
to our final objective, i.e. the coefficients A m and B m . The only non-vanishing coefficients 
are 

A m=0 = — {—2 lncosh(gar) — ln[l + i(b/a) tanh(gar)] — ln[l — i{b/a) tanh(gar)]} (37) 
nq 

B m =o = —(-l) q+1 {2iqbr - ln[l + i(6/o) tanh(gar)] + ln[l - i(b/a) tanh(gar)]} (38) 
7rq 

~ tanh J (qar) ( (1 — ib/aY 

m=qi = vq~3 Ul + z(-l)" +1 (6/a)tanh(gar)]J + 



[l + ib/a) j 



B m=q j 



[1 - i(-l)i +l {b/a) tanh(gar)]^ 
itanh J (gar) f (1 — ib/ay 



(39) 



nqj \[1 + i(-iy +1 (b/a) tanh(qar)y 

(l + ib/a) j 



[1 - i(-l)i +1 (b/a) tanh(gar)]^' 



(40) 



where j £ N. Note that since A m<0 = B m<0 = 0, the original coefficients A m and B m are 
equal to the A m 's and -B m 's respectively. 

The expressions for the RG-improved coefficients ensure that the solution for G and F 
does not depend on r. We still have the freedom to choose the arbitrary time r at will, and 
the obvious choice is t — r, since in this way the secular terms proportional to t — r dissapear. 
Given the RG-improved coefficients, we still have to plug them into Eq.(p7|) and perform 
the necessary summations to finally get the RG-improved solutions G(t,t) and F(t,t). 
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For a general dephasing, the resulting expressions are rather lengthy, so in the next two 
sections we will concentrate on particular cases. Firstly we study the case of dephasing 
= 0, which corresponds to translational modes, and dephasing = it, which corresponds 
to breathing modes. Secondly we analyze a case with only one mirror oscillating, similar to 
the one we studied in flIS| , but with a dephasing = 7r/2, which gives cualitatively different 
results. 



IV. TRANSLATIONAL AND BREATHING MODES 

In the present Section we consider that the cavity has translational modes (0 = 0), or that 
it has breathing modes (0 = tt). For the particular case of equal amplitudes Al = Ar, the 
former type of motion corresponds to the cavity oscillating as a whole, with its mechanical 
length kept constant (pictorically called an "electromagnetic shaker" |^0[), while in the latter 
type of motion the mirrors oscillate symmetrically with respect to the center of the cavity, 
the mechanical length changing periodically (an "antishaker" ) . Both for translational and 
breathing modes the expressions for the coefficients in Eqs. ( p7| - |40|) simplify considerably 
because 6 = 0, and the summations to get the functions G and F are straightforward. 
Setting t — r in Eq.(|27D we get the RG-improved solutions 



G(t, t) = 1 - -Im ln[l + C + (1 - C)e iq7Tt/A ] + eg^ (t, t) + e<? (2) (f , t) (41) 
A 7rq 

F(t, t) = G(t) + 24 S in(^/A) 1 + c2 + (1 _ 2 c C 2)cosM/A) , (42) 

where we have defined ( = exp[2gat]. The (now nonperiodic) RG-improved functions 

g {1) (t,t) and g {2) (t,t) are 

9 (.» (t , t) = 8in( ^ /A) __j^__ (43) 

9 <« (M ) = T (-ir^^ ri /A) 1 + ( , + [1 _^ (ji , ;tr (44) 

where, in the last formula, the upper sign corresponds to = and the lower sign to = ir. 
This solution for G and F is cualitatively similar to the one we obtained in the case for one 
oscillating wall with zero dephasing ||16|| . For the same reasons described in that reference, 
both RG-improved nonperiodic functions give negligible corrections to G and F in the long 
time limit (e _1 C t/A C e ~ 2 ). However, they are crucial for the solution to satisfy the 
correct boundary conditions at short times (t/A <^ e _1 ). 

The energy density inside the cavity is given by Eq.(|9|) in terms of derivatives of G and 
F. Since these expressions involve second and third derivatives of these functions, and since 
there is an initial discontinuity of the velocities of the mirrors, the energy density will develop 
5 function singularities that will be infinitely reflected back and forth between the mirrors. 
In what follows we will ignore these singularities 0. 



4 These singularities are of course artifacts of the sudden approximation we are using to describe 
the motion of the mirrors at t = 0. 
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The structure of the electromagnetic field within the cavity at long times strongly depends 
on the relation among amplitudes, frequencies and dephasings. If these are such that the 
coefficient a is equal to zero (remember that for the motions considered in this Section the 
other coefficient b is always null), then there is destructive interference. For equal amplitudes 
Al = Ar, this takes place for even q and dephasing <p = 0, or for odd q and dephasing 
<p = tt: all RG-improved coefficients A n and B n are null, and there is no motion-induced 
radiation enhancement whatsoever. If, on the contrary, a ^ 0, then we have constructive 
interference, which, for Al = Ar, is maximal for odd q and <fi = 0, or for even q and 
= 0. Radiation enhancement takes place: the electromagnetic shaker and antishaker have 
"explosive cocktails" at long times. In particular, for q > 2, the RG solutions G(t,t) and 
F(t,t) develop a staircase form. Within regions of t between odd multiples of A, there are 
a total of q jumps located at values of t for which the argument of the logarithm in Eq.fliTD 
vanishes, i.e. cos(g7rt/A) = ±1, where the upper sign corresponds to a > and the lower 
one to a < 0. In Figure 1 we show the form of the functions G and F for short times and 
in Figure 2 for long times. Note that in the long time limit they are practically the same. 
The energy density builds up a number of q travelling wave packets which become narrower 
as exp(— 2q\a\t) and higher as exp(4g|a|t), so that the total energy inside the cavity grows 
like exp(2g|a|t) at the expense of the energy pumped into the system to keep the mirrors 
moving as predetermined. In Figure 3 the profile of the energy density inside the cavity at 
a fixed time is depicted. We compare the case of the shaker with that of a single oscillating 
mirror. The difference in height and width of the peaks between these two situations are 
due to the fact that the parameter a for the shaker is twice that of the single mirror. This 
reflects how the cavity can enhance vacuum radiation. 

A rather different picture appears when one considers the q = 1 case, which corresponds 
to an oscillation frequency equal to the lowest eigenfrecuency of the cavity. In this case the 
energy does not grow exponentially, but oscillates around the static Casimir value. 

Now we turn to calculate the number of motion-induced photons inside the cavity. We 
assume that before the mirrors started to move the state of the field was vacuum, and that 
at time t = T, when both walls come once again to rest (], we define a new vacuum, in 
which there is an amount of real photons given by the Bogoliubov coefficients Eq. (|T0|) . To 
calculate these coefficients at a time t > T we need to know the form for the function G and 
F in the corresponding time intervals as appear in the integral expression Eq.([T0|). 

To this end let us discuss briefly how the RG-improved solutions G(t, t) and F(t, t) match 
the solutions to the problem of motionless walls for t > T. The non-local structure of Moore's 
equation implies that the solution for F(t) is the RG one F(t,t) up to t < T, and that for 
G(t) is the RG one G(t, t) up to t < T + A. Also, evaluating Moore equation for times t > T 
it follows that F(t) = G(t) for t > T. Finally, for t > T + A both Moore's equations can 
be combined to obtain the usual equation for a static cavity, so G(t) = F(t) = t/A + A(i), 
where A(t) is a 2A-periodic function that we must determine. If due care is taken of the 
boundary conditions at the moment when the walls stop, it is easy to see that this function 
can be written by periodizing the RG-improved functions G(t,t) and F(t,t) as follows 



5 For the motions L(t) and R(t) we are considering in this Section, this happens for times T such 
that T/A = 2k/q, (k £ N). 
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(45) 



and A(t) = A(t + 2A). 

Having now the form of the solutions for times after the stopping of the walls, we can 
calculate the Bogoliubov coefficients for late times t/A > T/A ^> e _1 in a manner similar to 
that of [13|] . Let us split the solution G(t) in Eq.([4l"l) in the form G{t) = G s (t)+Gnp(t), where 
the first part is G s (t) = t/A — 2/(7rg)Imln[l + C + (1 — Qe iqwt ^ A ] and the last part involves the 
RG-improved functions and g( 2 \ From a graph of these nonperiodic functions and of 
the difference between G(t) and F(t) one can see that, for long times, they are all bounded, 
much smaller than G s , and that they are effectively zero except in small time intervals which 
tend to zero as time increases. Therefore we can drop their contribution in the imaginary 
exponents of the integral representation of the Bogoliubov coefficient, and get 

1 f t/A+1 

f3 nm (t, T) = -a / — / dxexp{ — m[nG s {Ax) + mi]}. (46) 

2 V n J t/K _ x 

The function G s has a first term, linear in time, and a second one, that for late times 
becomes an oscillating function, its period being 2A/q and the amplitude of its oscillations 
being independent of e. Then the Bogoliubov coefficient can be rewritten as follows [] 

1 / I- 1 r2/q 

P nm (+ T) = - . / — e" i?r(n+m)(t/A_1) e ~ i7r (n+m)(2/q)k / ^ xe -in[(n+m)x+nf(x)}^ 

2 V " J I) 

v fc=0 JU 

with f(x) = — 2/(g7r)Imln[l + ( + (1 — () exp^iqnx)}, and ( = exp(2gaT). To go further and 
be able to perform the integral we make a piecewiese linear approximation for the function 
/, valid for late times. We concentrate on the case a < 0, for which £ — > at late times Q. 
From the graph of f(x) one can see that it can be approximated by 

forO < x < i - S 

f(x) = { -M 1 ~ ^ [ x ~ \) te\-&<*<\ + * , (48) 




for 1 _|_ § < x < 2 

q — — q 

where 5 = 2y/(/(qTr). With this approximate form the integrals become trivial, and after 
neglecting the integral over the middle interval which is proportional to 5, one can get a 
closed expression for the Bogoliubov coefficient, valid for m5 1. For the particular case 
q = 2 we get 



6 As we have anticipated, the dependence of the Bogoliubov coefficient on the time t > T is just 
a phase. 

7 The case a > gives similar results for the amount of created particles. The technical difference 
is that since £ — > oo for late times, the approximate function is different. 
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Next we need to calculate the sum over n in order to find the amount of motion-induced 



photons in the n-th mode after the stopping time T. Using the summation formulas of JT3 

we get 



N m {T) ' 



m? 



.25) V2^' ■ <50) 

Recalling that 5 is a function of T, and taking the T-derivative we find the rate of photo- 
production 

^-^-(-m pi) 

These results are valid for eT/A ^> 1 and not for very large wave numbers mS <C 1 f| The 
number of photons per mode grows linearly in the stopping time and the rate, for late times, 
approaches an asymptotic value that depends on the value of a, i.e. one the relation among 
amplitudes, frequency and dephasing. Both for the shaker and the antishaker, motion- 
induced radiation is enhanced in comparison to the case of a single oscillating mirror in a 
cavity. Indeed, in the former cases the rate of photoproduction is twice that of the latter 
case. Photons are created in the odd modes only, whereas their amount in even mode is zero 
(it may be different from zero in the next-to- leading order approximation). This situation 
is typical of processes involving parametric excitations 13. 



V. A DEPHASED OSCILLATING BOUNDARY 



In this Section we discuss another particular motion of the walls, namely one for which 
the left mirror is static and the right one oscillates resonantly, with a dephasing <p = tt/2. The 
motion we consider in then R(t) = A — 2eAn sin 2 (qnt/ 2 A), which for q = 2 corresponds to 
the small e-expansion of an exact solution to Moore's equation studied in ||. Our motivation 
for studying this peculiar case is that for this motion we have a = and b ^ 0, which, as 
we have anticipated, gives cualitatively different physical results. 

The expressions for the functions G(t,t) and F(t,t) are obtained taking the limit a — > 
of Eqs.(ETJ|4y). The result is 

iqbte iqnt/A 



G(t,t) 
9 {1 \t,t) 



F(t,t) 


A R 



t 



-Im In 



7iq 



1 - 



1 + 



- i(-l)i +1 qbt 
l) q+l cos(qnt/A) 



+ eg®(t,t) 



A l + 2qbtsm(qirt/A) + 2(qbt) 2 [l + (-1) 9+1 cos(gvrt/A)] 



(52) 
(53) 
(54) 



8 There is a further restriction that comes from the fact that we are using a sudden approximation 
for the motion of the mirrors at t = and t = T. Indeed, if we assume that t s is the characteristic 
time for the mirror to come to rest, the sudden approximation will be valid for modes such that 
m < A/r s . 
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where A e g = A(l — sAr/A) is the time-averaged length of the cavity for t > 0. 

The solution G(t, t) develops a staircase profile, the jumps being located at values of t for 
which the argument of the logarithm in Eq.(p2[) vanishes, i.e. for cos(g7rt/A) = ±1, where 
the plus sign corresponds to even q and the minus sign to odd q. The energy density for 
this type of motions also consists of a series of q peaks that travel between the mirrors. The 
cualitative difference is that in this case the height of the peaks grows as (qbt) A , their width 
decreases as (qbt)~ 2 and the total energy contained in the cavity grows quadratically rather 
than exponentially. This follows from the fact that time enters in the logarithm of Eq. (p2[) 
as a power law instead of an exponential, as in Eqs.(|I],I2). 



Next we calculate the amount of mot ion- induced radiation for this case. To this end 
we assume that at time t = T the wall comes to rest, R(t) = A e g for t > T, where T is 
of the form T/A = (2k + l)/(2q), (k e N). This choice for the stopping time simplifies 
the computation of the the Bogoliubov coefficients (3 nm . In such a case Eq.(f46l) is slightly 
modified 

1 [m f t/A cs +1 

[3 nm (t,T) = — dxexp{-i7r[nG s (A ea x) +mx]}. (55) 

1 V n Jt/A cS -i 

Now G s consists of the first two terms of Eq . (|52D , the first being linear in time and the second 
one being an oscillating function for late times, whose period is 2A/q. The Bogoliubov 
coefficient can also be rewritten in a way similar to Eq.( 



1 I 5~- r2/q 

(3 nm (t, T) = -J — e -i*(n+™)W^&-l) e -Mn+ m )(2/ q )k / ^ e ^[{n+m)x+nf{x)] /ggN 

2Vn ^ J Q 

since in the interval [t/A eS — 1, t/A eff + 1] there are a total of q(A cS / A) rs q periods. Here 

iqbA eS xe iq7TxA ^/ A 



2 

fix) = Imln , 

JK ' nq [ l-i{-l) c '+ 1 qbA cS x_ 
The piecewiese linear approximation for the function / is in this case 

^ - x - i forO < x < 5 

2 I q — — 



(57) 




qS) for5 <x<l 



where 5 = [g7r(g6T) 2 ] _1 ^ 1 for late times. Now the integral in Eq.(p6[) is straightforward, 
and after dropping the integral over the first interval [0,5], which is proportional to 5, one 
can get a closed expression for the Bogoliubov coefficients, valid as long as m5 1. For 
comparison with M we concentrate on the case q — 2. In this case we have 



|2 _ 2m. . m+nl l-(-ircos(37rn5/2) 



\(3 nm (T)\ 2 = —[1 + (-l) m+n ] ' ' (59) 

rm z (m + dimo / 2) z 



'We concentrate on even frequencies, for which b < 0. For odd frequencies, the results are similar. 
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Finally, we perform the summation over n to get the number of created photons in the n-th 
mode. Using the same summation formulas as in [13| we get 



N m (T) 



' (2m 



;-l) m ln 



3tt5 



(60) 



Replacing the value for 6 and taking the T-derivative, we get the following formula for the 
rate of photon production, valid in the limits m5 <C 1 and eT/A ^> 1 



dN m {T) 
dT 



rriTC- 



'[I 



(61) 



We see that the number of photon per mode grows logarithmically in the stopping time, and 
in consequence the rate of photon creation decreases towards zero. Similarly to the case of 
the vibrating cavity, photons are produced only in odd modes. 



VI. CONCLUSIONS 

In this paper we have presented a unified and analytic treatment of the dynamical Casimir 
effect in a one dimensional resonantly oscillating cavity for arbitrary amplitudes and de- 
phasings. We have derived a generalization of Moore's equation to describe the state of the 
electromagnetic field inside the cavity with two moving mirrors. Using a technique inspired 
in the renormalization group method, we have found a solution to the set of generalized 
Moore's equations which is valid both for short and long times. The physical behaviour 
of the moving cavity depends crucially on the relation among amplitudes, frequency and 
dephasing. We have shown that for certain cases there is destructive interference and no 
radiation is generated. For others, there is constructive interference and motion-induced 
photons appear. When this takes place, the way the energy within the cavity and the num- 
ber of created photons grow in time depends on the relation among the above variables. For 
certain motions the growth of the energy density is exponentially and for some others it is 
a power law. 

We hope in the future to apply the RG method to more realistic situations, such as three 
dimensional oscillating cavities with rectangular or spherical shapes. 
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Note added: While we were finishing this manuscript we received a paper |21| where the 
problem of photon creation in a cavity with two moving mirrors is analyzed using a different 
method. It is shown that if the frequency of the vibrations is not exactly a resonant one, 
photoproduction is highly supressed for strong detuning. 
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FIGURES 




FIG. 1. G(t) and F(t) vs. t/A as given by Eqs.(||j4|) for small times q\a\t/A < 1. Note that 
the function F(t) departs from the straight line at t = 0, while the function G(t) departs from it 
at t = A, as dictated by the initial boundary conditions. The parameters are those for a shaker 
with A L /A = A R /A = 1, q = 3, <j) = and e = 0.03. 
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FIG. 2. G(t) and F(t) vs. t/A as given by Eqs.@g|) for long times q\a\t/A > 1. At these 
times both functions coincide and take a staircase profile. The parameters are the same as in figure 
1. 
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FIG. 3. Energy density profile between plates for fixed time t/A = 15.4. The solid line cor- 
responds to a shaker with Al/A = An/ A = 1, q = 3, <j) = and e = 0.01. The dashed line 
corresponds to a cavity with a single oscillating mirror {Al = 0) with the same parameters. Note 
that the height of the peaks for the shaker grow as exp(47rqet/A), while for the single mirror as 
exp(27rgei/A). 
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